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Abstract 

The properties of artificial isotropic microwave composites which would possess simultane- 
ously negative permittivity and permeability are studied theoretically. Four kinds of com- 
posites are considered. Two of them concern media with split-ring resonator particles with 
different particle arrangements. The other two are realized with Omega particles. An ana- 
lytical antenna model of the electromagnetic behavior of a split-ring resonator (SRR) is sug- 
gested and verified by numerical simulations. Next, material parameters of composite media 
made with SRR or Omega particles are calculated. It is shown that the Omega-composites 
are more prospective for obtaining negative real values of permittivity and permeability than 
the split-ring resonators. 

1 Introduction 

V. G. Veselago studied in sixties [1] the electromagnetic properties of materials whose per- 
mittivity and permeability have negative real parts simultaneously. He showed that these 
media exhibit unusual effects like anomalous negative refraction of electromagnetic waves and 
existence of backward waves (when the group velocity and the Poynting vector are directed 
opposite with respect to the phase velocity). This concept implies many other interesting 
features of wave propagation that were reviewed in [1]. 

These Veselago media do not exist in nature. The interest in creating such media artifi- 
cially, for optical applications, has been indicated by J. Pendry in [2]. This kind of medium 
can be obtained with arrays of conducting inclusions of a special shape. It is obvious that 
the response of these materials will be resonant. Therefore, negative refraction is possible 
only within a more or less narrow frequency band (or bands) [3]. This "negative" band can 
be a part of a resonant band of inclusions that are used in the composite. This is quite clear 
analyzing the Lorentz model of the resonant dispersion. 

At this stage of studies, since it is easier to use resonant inclusions with millimeter 
dimensions, it seems to be reasonable to focus on applications for the microwave region of 
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the spectrum (optic composites will need micrometer particles) and to study their response 
experimentally. How to prepare these composites in a best way? Which shape of inclusions 
is more suitable? These are the general questions addressed by this article. 

Works that have theoretically predicted negative permittivity for composite media have 
been reviewed in [4]. However, publications containing the same results for the permeability 
are not so common. In [5] and [6] it was indicated that it was a very difficult problem to 
obtain a magnetic resonant response of inclusions (in order to get a negative real part of 
/i e //). Probably it was the reason for which these authors have chosen a rather complicated 
particle named as split-ring resonator (SRR), see Fig. 1. They maintained in these works, 
that SRR has only a resonant magnetic polarizability and no resonant electric one. 

In [6] it was suggested that to obtain negative permittivity simultaneously with negative 
permeability. SRR inclusions must be combined with a lattice of straight wires. This complex 
system is a very cumbersome technical solution. The only special case where wave propagates 
normally to the wires (which corresponds to [7]) is suitable to observe the negative refraction 
in such a structure. It is clear that the structures from [6] cannot be considered as an isotropic 
backward-wave composite. 

In this paper we show that the SRR particles alone allow to obtain simultaneously neg- 
ative permeability and permittivity. Indeed, these inclusions also possess a resonant electric 
polarizability that has not been discovered in [5] and [6]. We also point out that the analyt- 
ical model of SRR suggested in [5] does not give a satisfactory explanation of the frequency 
behavior of effective permeability. Moreover, it is possible to realize an isotropic composite 
with negative refraction choosing properly the geometry, arrangement and concentration of 
SRR particles in a host dielectric matrix. 

However, we have found that SRR is not the best choice for obtaining a backward-wave 
medium at microwaves. Actually, the resonant electric response of SRR is weaker than 
the magnetic one. We show that it is better to use an inclusion known as Omega particle 
(or fi-shaped conducting particle) suggested in 1992 by Saadoum and Engheta [8] since its 
magnetic and electric responses of are of the same order. 

Numerical and analytical modeling of Omega particles have been recently made in [9]- 
[12] and the properties of uniaxial and planar Omega composite media have been derived. 
It appears that these media are prospective for creating low-reflecting shields and antenna 
radomes in microwave applications. With this kind of objectives, in 1995-1996, we studied 
theoretically uniaxial Omega composites and noticed that sometimes the real parts of both 
permittivity and permeability became negative, especially in the cases with high densities 
of particles. Thus, as it was not suitable for our goals, the experimental samples of Omega 
composites which were prepared under guidance of S.A. Tretyakov [9] were dilute mixtures in 
order to avoid the apparition of negative properties. Negative permittivity and permeability 
of Omega composites has been predicted in the paper [13]. However, the model of the 
material parameters presented in this work and based on the transmission-line model of 
Omega particle was contradictory (particle size was taken of the order A but the constitutive 
relations were introduced as if the composites were continuous media) and violated some 
physical foundations (Kramers-Kronig relations, etc). 

In the following we consider the effective properties of isotropic composites of both Omega 
and SRR particles embedded in a dielectric host medium. We study two different kinds of 
particles arrangement. The first one is a mixture with a random orientations of particles. 
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The second one is obtained by positioning the particles at the sides of a cubic unit cell, like 
it is shown in Figs 2 and 3. Different inclusions and structures are finally compared. 

2 Theory 

Both SRR particles and Omega particles are bianisotropic particles. In other words, an 
external electric field produces a magnetic dipole in the particle and, at the same time, an 
external magnetic field produces an electric dipole. For SRR, this magnetoelectric effect is 
reduced by choosing the opposite positioning of the splits of both loops. However, it cannot 
cancel out exactly since the loops are of different sizes. 

Let a SRR particle be located in an infinite and homogeneous dielectric medium with 
permittivity e. Consider the electromagnetic response of an individual SRR particle to local 
magnetic and electric fields. The polarizabilities of a SRR particle are dyadic coefficients re- 
lating the local electric and magnetic fields E loc , H loc with the induced electric and magnetic 
dipole moments: 

p = ^ e -E loc + ^ em -H loc (I) 
m = ^ me -E loc +^ mm -H loc (2) 

The electric and magnetic moments of the SRR as a whole are vector sum of the moments 
induced in the outer and inner loops (indices 1 and 2, respectively) : 

P = Pi + P2 (3) 
m = in, |- m 2 (4) 

In this section we evaluate analytically the polarizabilities of a SRR particle, assuming that 
the inclusion is fabricated from a round wire of radius r$. The section of the wire is small 
enough to satisfy the inequality r <C ai,2, where a ii2 are the averaged radiuses of rings 1 
and 2 (see Fig. fa). Distance 5 between the broken ends of each ring is assumed to be small 
compared to a 12 . The analytical model of the Omega particle will not be reproduced here. 
Analytical antenna theory and its verification can be found in [10]. 

Finally, using the calculated polarizabilities we derive for each kind of particles the ma- 
terial parameters through the well-known Maxwell Garnett model. 

2.1 SRR particle eigenfrequencies 

We consider the SRR particle as a pair of two mutually coupled resonant scatterers. The 
complex amplitudes of induced currents at the splits of both rings can be written as follows: 

h(Z cl + Z L1 ) + I 2 Z 12 = £ H i + £ei, (5) 
h(Zc2 + Z L2 ) + I\Z\2 = £h2 + £e2 (6) 

Indices 1 and 2 refer to the outer and inner ring respectively. £ H and £e are the electromotive 
forces induced by the external fields. I± j2 are the currents at the splits of the rings. Z L ^ L2 
represent the inductive loop impedances (Zli,l2 = j^-^1,2), and Zci,C2 are the capacitive 
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impedance {Zci,c2 = l/(j^Ci,2)) of the slits. For a magnetic field excitation we can neglect 
the individual loop capacitance which has only a small influence. For the study of electric 
polarization, influence of the proper capacitance of one loop must be taken into account. 
For that we use the known theory of a split loop antenna [16] and [15] which takes into 
account the capacitive properties of an individual loop scatterer and its radiation losses. 
Z12 in (5), (6) is the mutual impedance of broken loops. We have chosen to include directly 
the capacitive mutual coupling in the term C\^- Otherwise the analytical model of mutual 
coupling would become too involved. With these considerations the mutual impedance Z\ 2 
in (5), (6) is expressed only as a mutual inductance M (Z 12 = jivM). 

The capacitances C\^. are composed of two contributions. The first part C$ originates 
from the splits of each loop. The second one is due to capacitive coupling of loops and can 
be seen as a consequence of the following effect. Consider first as a reference the conductor 
of one loop (the outer loop of a SRR particle, for example). Next, place another conductor 
parallel to the first one (in our case, the inner loop, with the split at the opposite side). Then, 
assume that one half of the first loop (e.g. the upper part) contains positive charge +Q, and 
the second half of the first loop (e.g. the lower part) contains negative charge —Q. These 
charges are distributed somehow along the loop perimeter. As the two loops are positioned 
closely, in the upper half of the second loop, a negative charge — q is induced. Hence, the 
distribution of the global charge — q around the second loop repeats the distribution of +Q 
in the first loop. Note that we do not have necessarily the correspondence Q = q. Of course, 
the same thing happens in the lower part of the structure, where the induced charge +q 
appears. So two equivalent capacitances C mut arise from this situation as shown in Fig. 4. 
Hence, we obtain an effective additional capacitance which is equal to C mut /2, connected in 
parallel to the proper capacitance of the split. One can note that this equivalent scheme can 
be improved taking into account that there is an additional capacitance of the second slit. 
However, we are trying to derive a model of SRR as simple as possible and the complexity, 
increased by this additional capacitance, is not justified. 

Values of inductances L 12 , mutual inductance M and total capacitances Ci j2 are analyt- 
ically calculated in Appendix. 

To obtain the eigenfrequencies of the SRR particle, we let the right-hand side of (5), (6) 
vanish and equate the determinant to zero: 



where r 12 are small resistances that represent the radiation and the Joule losses of both 
loops. They lead to a small imaginary correction term denoted as k{uj) in equation (9): 



If we assume that the SRR particle is perfectly conducting, the term k describes only the 
radiation losses of a magnetic dipole and it can be derived from the energy balance condition, 



(Zci + Z L i)(Z C 2 + Z L2 i) — Z\ 2 = 



(7) 



Equation (7) can be written as 





{u 2 -u 2 _){u 2 -u 2 + ) +jk(u) = 



(9) 
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as it was done in [17] for arbitrary dipole scatterers. Comparing (8) with (9) we find the 
eigenfrequencies: 

'uj±\ 2 uxL 



+ 



uj 2 



2 



\ 



-4 



(10) 



where we denote 



1 



^C X C 2 {L X L 2 -M 2 ) 



1,2 



Cl,2-^1,2 



Frequencies (10) correspond to two eigenmodes. If the frequencies c<j_ and u + differ 
dramatically the two eigenmodes of the induced current are excited at different frequencies. 
This leads to two separate resonances (that magnetic and that electric ones) of the particle. 
However, we will see below that in practical cases u + ~ c<j_. The two resonance bands 
overlap, and we obtain a result close to that of the Lorentz model with only one resonance, 
as it goes for simple scatterers. As a consequence, it matches qualitatively the results of [5] . 



2.2 Magnetic polarizability of SRR particle 

Next, we are going to derive the polarizabilities of SRR particle using (5), (6) and definitions 
(1) and (2). First, we study the magnetic polarizability. For a magnetic excitation (E loc = 
0, H loc = z H loc , where the z-axis is perpendicular to the loop plane as it is shown in Fig. 
2) we obtain: 

£m,2 = -jufioS lt2 H loc , m lj2 = h,2l^oSi, 2 (11) 

Here Si^ = na\ 2 are the effective loop areas. 

The only component of the magnetic moment is vertical, i.e. 

a mm = z O z a mra 



Then, relations (11) together with (l)-(6) and (7)-(9) lead to the following expression: 

A+jK 



a mm - A , 77. \ lZ ) 



where we have denoted 

A = {u 2 - u 2 _){uu 2 - u 2 + ) (13) 
B 1 = S 1 (S 1 L 2 - MS 2 ), B 2 = S^S^ - MS 1 ) 



D- & 



L X L 2 - M 2 
and 

= (l %rj (14) 

Now let us find the term k. The theory from [17] based on the balance of received and re- 
radiated power for any particle with magnetic dipole polarizability a^ m gives the equation 

1 P 

Im— = -A- (15) 
C 6vr/i 
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Comparing this formula with (12) we come to the following expression for the term k: 



u?Dk* [B^u 2 - u> 2 ) + B 2 (u 2 - uj> 2 )} 
K = (16) 

67T/i 

With an appropriate calculation of inductances and capacitances (see Appendix), it is now 
possible to completely evaluate the magnetic polarizability a z r 



^zz 



2.3 Electric polarizability of SRR particle 

Now we are interested in the electric polarizability. Consider the action of local electric fields 
on a SRR. Two possibilities exist to induce polarization on the SRR particle. 

First, we study the case where the electric field is along the horizontal direction (x— direction 
for an individual particle). The only electro-dipole mode of the current is excited in both 
loops in this case. The induced charges are concentrated at the splits and at the oppo- 
site points and the polarizability component a** is quasi-static. It does not depend on the 
presence or absence of the split and equations (5), (6) are useless in this case 1 . 

In [14] we derived formula (28) for xx— component of the polarizability of a circular loop 
with area S = ira 2 . Denoting this polarizability as ciqs we obtain from [14]: 

ASJ'^ka) 
<: " oor]Ai{ka) 



In this equation, rj = ypio/tto is the wave impedance of the host medium, k = Uy/fi eeo is 
the wave number in the medium, J[(ka) is the derivative of the Bessel function of the first 
order, Ai(ka) is one of so-called King's coefficient known from the theory of loop antennas 
[15], [16]. King's coefficients with a very high accuracy are approximated in [15] and take the 
following form: 



Ai(ka) 




O.QQ7(ka) 3 - 0.207(ka) 1 



IT 

-j0.333(ka) 2 - j(0.133ka) 4 + j0.026(kaf 
Mka) = £ (l 0g !2 _ 2 ) + «>• W - °- 267 <fa) 5 _ m67[kar _ j0 .033( M < 

Mka) = ( k ^A) ( log *L _ 0^667) + ~ 0Aka + °- 21 ' fa » 3 " ° ° 86 < fa » 5 



7T / \ r n / 7T 

\4 -n m on,.„\6 



-j0.05(A;a) 4 - j0.012(fca) t 
Since we have two loops, we must sum up their ciqs to obtain the result: 

aXX = 4S 1 J((fco 1 ) _ AS 2 J[(ka 2 ) 
ee " urjAiikai) ujr]A 1 (ka 2 ) 



1 Because — and the electromotive force £e induced at the split is also equal to zero: charges are 
the same at the edges of each split. 
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Now we consider the case where the electric field is directed along y— axis. Here the resonant 
electric polarization appears in both loops of the SRR particle. To obtain the resulting 
polarizability, we combine the theory of a single loop with resonant electric polarization (see 
[16] and [14]) with equations (5) and (6). 

In [14] the following formula has been derived for the yy— component of the electric 
polarizability of a single loaded loop of radius a: 

= A„a>J>(ka) ( _2j l\ 

Here Y sp u t is the admittance that is connected to the loop split. For both loops we have 
^spiit = V-^spiit = 3^^1,2 which include the mutual capacitance of the loops. Y L = l/Z L 
is the proper loop admittance. In our case of coupled loops, we should replace Y L by the 
value Y"i j2 which takes into account their mutual inductance. For loop 1, we can say that it is 
composed of the proper impedance of the loop Zli and of the additional impedance Zuh/h, 
which is the contribution of loop 2 into loop 1. Then, for loop 1, we have to substitute into 
(18): 

— = Z\ = Zli + Z12 — 
i\ *i 

Similarly, for the second loop we have: 

1 J i 

— = Z2 = Zl2 + Z12 — 
12 12 

Denoting the ratio /1//2 as £ and taking into account that Z 12 = jujM, we come to the 
following formulas for the yy— polarizabilities of both loops: 

a vv _ 4^j(fcai) L 2j 1 \ 

el urjA^kai) \ 7rrjA 1 (ka 1 )Y 1 + juCiJ 

4^2/ 2J 1 \ 

oo7 1 A l (ka2) V 7Tr]A l {ka 2 )Y 2 + juC 2 J 



where 

Yx 



Y, 



Ll 



1 , jujMY L1 



1 + jujMY L2 Z 
The proper loop admittance is given by [15]: 



1/1 2 2 \ 

L ~J^I \A (ka) + A^ka) + A 2 {ka) ) ^ 



Expressions of Y L1L2 follow from (21) with substitutions a = a 12 . 
Equations (5), (6) allow to find £: 

^ _ £ei(Zl2 + Zq2) — £E2Z\2 
£e2(Zl! + Zci) — £e\Z\2 



(22) 
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If a loop with proper impedance loaded with the impedance of the split Z sp u t is submitted 
to an external electric field, an electromotive force is induced. Its value is given by formula 
(26) from [14]: 

£e ~ piMka) [ 1 + ^nT^) E « (23) 

Substituting the values a\^, and r into (21) and (23) allows to calculate the electromotive 
forces induced by an external electric field in each loop of the SRR particle. Then, substitut- 
ing expressions of £ei,e2 into (22) with the help of relation (21) we can find £. Consequently, 
we derive an expression of Y 12 which is then substituted into (19) and (20). Finally, we 
evaluate a y e y e as the sum 

all = all + all (24) 

At this point, we have at our disposal a complete analytical model of the electric and 
magnetic polarizabilities of an individual SRR particle (other components of a ee and a 
are zeros). 



mm 



2.4 Effective polarizabilities of a planar bianisotropic particle in 
isotropic mixtures 

An isotropic mixture with a random orientation of particles can be considered as a mixture 
realised with isotropic particles which have averaged scalar polarizabilities [19]: 

av _ "ee av _ a mm av _ av _ n /r>,-\ 

"ee — q ' a mm ~ q ' V ~ "me ~ U l Z0 J 

a=l ° a=l ° 

Here we have designated Cartesian components with index a — x,y, z. The particles under 
consideration do not have cross components: a^f = for a ^ (3. Moreover, we have seen 
with our analytical model, that the only polarizabilities that are non-zero are: , a vy e and 

"'mm' 

If SRR or Omega particles are arranged on the sides of cubic unit cells (see Figs. 2 and 
3), the magnetoelectric effect also disappears and dyadics a em) a me of the individual particles 
do not contribute into material parameters. Furthermore, mutual coupling of six particles 
belonging to the same cell is weak enough if the cell size d ce u is a bit more than the particle 
size dp (the particles do not touch one another). Now, let us evaluate the equivalent effective 
polarizability of each kind of unit cell ignoring the mutual coupling between the particles of 
a cubic cell. 



2.4.1 Effective polarizabilities of a cubic cell of SRR particles 

The cubic cell under consideration is presented Fig. 2. An electric field directed along an 
arbitrary edge of a cubic cell excites the electric dipoles in four SRR particles. For two 
particles E loc is parallel to the line connecting the rings splits and corresponds to the polar- 
izability cigg. For the other two particles E loc is perpendicular to this line and corresponds 

to a yj. 



8 



Consider now the magnetic field effect on the cell. A magnetic field directed along an 
arbitrary edge of a cell induces two equivalent magnetic moments in two opposite particles. 
So, the magnetic polarizability of a unit cell turns out to be the double of af^ m . 

To summarize, the cubic unit cell is equivalent to an isotropic particle with electric and 
magnetic polarizabilities given by 

a' ee = 2a~ + 2a™ a'mm = 2<C W 

2.4.2 Effective polarizabilities of a cubic cell of Omega particles 

For a unit cubic cell of Omega particles (see Fig. 3) the situation is quite different. Here 
the resonant excitation by the electric field directed along the arms of Q (along y-axis in 
the model of an individual particle) is mainly due to the presence of the arms [10]. There- 
fore the total electric polarization of two opposite omega particles keeps resonant behavior. 
Within the resonant band the quasi-static component of Omega particle polarizability a** n 
describing the response of a particle to the electric field applied normally to the arms (see 
Fig. 1) is small compared to a v e v e n . Then one obtains for a unit cell from Omega particles 
approximate relations 

a' ee « 2<f, a' mm = 2<C (27) 
where a^ n , can be taken from [10]. 

2.4.3 Material parameters of composites 

The material parameters of the isotropic mixture are given by the standard Maxwell Garnett 
formulas [19]: 

'No' N 2 o' a' \ 

(28) 



~eff 



ee 



In these equations 




3eo/io 



N 2 a' a' \ 

H eff = 1 + F ( 3e)i e ) ^ 

Y'-V n' 

^i — l ee mm _|_ "ee^mm 

3eoe 3/io 9eoHQe 

These equations are general for the different kinds of composites we are studying. De- 
pending on the composition of mixture, only expressions for polarizabilities change. 



for random media, which are realized by mixing particles in an isotropic and random 
way in a host medium, N denotes the concentration of particles, and a' ee = a™, a' mm = 
a mm are S iven b y ( 25 ) 



for cubic-cell media, which are obtained with previously mentioned cubic unit cells, iV 
denotes the concentration of cells, and a' ee and a' mm are given by (26) or (27). 
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3 Model validation for SRR particle 



Our analytical model is below compared with numerical simulations. Numerical model which 
allows to simulate the averaged polarizabilities of an arbitrary bianisotropic particle with thin 
wire of round cross section, that allows to calculate also material parameters of a random 
composite has been derived from [20], [21]. Since its publication this method has been widely 
used to simulate bianisotropic composites. 

We compared the analytical and numerical models of the SRR prepared with wires of 
radius r = 0.2 mm. Numerical example presented here corresponds to the following param- 
eters. The averaged radiuses of the inner and outer loops are equal to a 2 = 2.1 and a\ = 2.7 
mm. The length of the broken part is taken 5 = 0.1 mm or 5 = 0.4 mm (5 is the same in 
both loops). The host medium permittivity e = 1.5. 

In Fig. 5 the real part of a^ m is shown calculated in both analytical and numerical ways 
for the case 5 = 0.1 mm. The imaginary parts of a™ m for this case are shown in Fig. 6. 
The real and imaginary parts of af e calculated in both ways are shown in Figs. 7 and 8, 
respectively. 

Comparing the resonant frequencies of a™ and a™ m , we can note that the difference is 
quite negligible. This means that the resonant bands overlap (u + c<j_) as we mentioned 
previously. 

We can observe also, that the analytical model shows a weak additional resonance for 
the electric polarization around 5.2 GHz. This is a defect of our analytical theory which is 
related to the adopted approximations. We think that the second resonance follows from 
our model of loops mutual coupling which is not accurate enough. Adding a very small 
imaginary value to the right-hand side of formula (31) we completely suppress this second 
resonance without changing the values of a mm and a ee within the first resonant band. Of 
course, this correction has nothing to do with the theory but it confirms that the model of 
Zi t 2 is not accurate enough and can be improved. 

The electric polarizability calculated in both analytical and numerical ways for the case 
5 = 0.4 mm (the other parameters are the same as above) is presented in Fig. 9. 

The agreement between two models is quite the same for 5 = 0.1 and 5 = 0.4 mm. The 
same remarks could be made for a mm . 

Polarizabilities of SRR particle are very sensitive to deviations of the particle parameters, 
and we can consider that the agreement between the two models is satisfactory. It makes 
our analytical theory suitable to calculate the effective composite parameters and, hence, 
allows to solve explicitly the inverse problem for such composites. 

4 SRR and Omega isotropic composites 

In this section we analyze the material parameters of effective media obtained with SRR 
particles or Omega particles. The mixture is fabricated by 

• distributing the particles with random orientations in a host medium 

• positioning the particles at the sides of a cubic unit cell (filled by the same medium as 
the matrix) as shown in Fig. 2 and 3. 
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4.1 SRR composites 



In Figs. 10 and 11 both e e // and /i e // of a random SRR composite are shown (calculated 
in both analytical and numerical ways). Particles geometry and background permittivity 
has been taken the same as above (the split width is 5 = 0.4 mm). The concentration of 
particles has been taken N — 1 particle/cm 3 . For the particle maximal size 5.8 mm this 
concentration allows to use the Maxwell Garnett averaging procedure. 

We can observe on Fig. 10 that the permittivity is not becoming negative within the reso- 
nant band (whereas both models of the permeability allow the negative values for Re(/z e //)). 
Here, the permittivity of SRR composite remains positive because of the restriction of dilute 
mixtures (imposed by the validity of the Maxwell Garnett model). 

In the case of SRR cubic media, since we apply the Maxwell Garnett model to the unit 
cells considering them as particles, we can take the concentration of SRR particles greater 
than in the case of a random mixture, and that, without violating the restriction of dilute 
mixtures. 

As a consequence, in Fig. 12 we can see that Re(/z e //) and Re(e e //) are negative simulta- 
neously. The concentration of cells has been calculated for the case of unit cells arrangement 
as presented in Fig. 3. The size of the cell is equal to 6 x 6 x 6 mm 3 , and the averaged 
distance between the centers of two adjacent cells is close to 11 mm. This allows to consider 
the composite with such cells as a dilute mixture. 

4.2 Omega composites 

In Fig. 13 and 14 we present the results calculated within the frame of the Maxwell Garnett 
model for an Omega composite prepared from cubic units as shown in Fig. 3. In that 
case, the maximal size of the Omega particle is equal to 6.4 mm. The cell size is equal to 
6.5 x 6.5 x 6.5 mm 3 . For an illustration, we have chosen the following dimensions. The 
radius of the loop is 3 mm and the radius of the wire is r = 0.1 mm. The length of an arm 
is L = 3 mm. The host medium permittivity is e = 2. Note that for cases when 5 < a, 
parameter 5 does not influence the polarizabilities at all (see [10]). 

In Fig. 13, the real parts of the permittivity and permeability of the Omega composite 
versus the frequency are shown. In Fig. 14, the imaginary parts of permittivity and per- 
meability are presented. We can see that Re(e e //) and R,e(fi e ff) reach —1 practically at the 
same frequency. This is a very exciting result if we refer to paper [2] where the perfect lens 
has been predicted for the case i?e(e e //) = Re(fi e ff) = —1. 

The case of random Omega composite (Fig. 15) with density N = 1 particle/cm 3 only 
differs by the heights of the resonant peaks of the permittivity. 

5 Conclusion 

Based on the analysis of electromagnetic phenomena that appear in a split ring resonator, 
we have developed an analytical model of SRR particles. We have shown that the SRR 
particles can exhibit, on their own, simultaneously magnetic and electric polarizability. The 
proposed model has been successfully compared with a numerical approach. So, these results 
point out that previous works on SRR do not give a satisfactory explanation of the behavior 
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of this particle. Next, we have calculated with the Maxwell Garnett approach the effective 
parameters of composites realized with SRR particles. We have shown that it is possible 
to obtain an isotropic composite with negative refraction choosing properly the geometry, 
arrangement and concentration of SRR particles in a host dielectric matrix. The way of 
making these composites has been presented. With this technique, we have obtained an 
effective material with simultaneously negative real part of permittivity and permeability. 

However, we have shown that SRR is not the only and not the best choice for obtaining 
a backward-wave medium at microwaves. Actually, we can obtain the same kind of medium, 
using Omega particle instead of SRR particles. With these inclusions we have derived a 
material which gives Re(e e jj) = Re(/z e /j) = —1 in the resonant band. This is the wished 
result exposed by the inventor of the SRR particle. Finally, we can conclude that the Omega 
particle is a more appropriate choice in order to create microwave composites with negative 
refraction. 

We have to point out yet, that our model exhibits a high level of radiation losses within 
the resonant band. As a consequence, if material samples are so lossy, it would be problematic 
to measure the negative refraction on it. These high radiation losses are due to the Maxwell 
Garnett model which implies a random organization of particles and is consistent with the 
Kramers Kronig relations. However, the proposed cubic composite can be classified in a 
middle position between random media and lattices. The composite is constructed as a 
lattice, stacking up cubes in a tree dimensional chessboard assembly, but it is not a perfect 
lattice. There is indeed, a certain random organization in the positions of the particles on 
the sides of the cubic cells (as it can be seen on figures 2 and 3). So, we can reasonably 
expect that losses given by the Maxwell Garnett model will be overestimated. 
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6 Appendix 

Inductance of a round loop of radius R made with a wire of cross section radius r can be 
found e.g. in [16]: 

L = MR (log (^) - 2) (30) 

Inductances of both rings of SRR can be obtained with (30) with the substitution R = ai >2 . 

Formula (5.23) from [18] for two concentric loops whose radii a±,a 2 are close to one 
another, gives the following approximation for the mutual inductance: 



M = fx ai 



x)logQ) -2 + C 



(31) 



where x = ( a i ~ a 2 )/2ai is a small parameter (0 < x < 0.25). In (31) we have neglected the 
terms of order £ 2 . 

Now, consider the capacitances included in the loops splits. We can write in accordance 
with Fig. 4 

C l = C sl + ^f, C 2 = C s2 + ^f (32) 
For C s we have the simple equation of a conventional (parallel-plate) capacitor 

C = (33) 

if the split width 5 is much smaller than the wire cross section diameter. Otherwise, if 5 ^> ro 
one has C s i ;2 C mut and the capacitance of the split can be ignored. We do not study the 
cases when r ~ S. 

Next, we calculate C mut . The capacitance per unit length of two parallel wires of diameter 
2r with the distance between the wire axes d = a 2 — a\ denoted below as Co is known [22]: 

arccosh(^) ( 

Roughly (but successfully) we can consider C mu t as the capacitance between two equivalent 
parallel straight wires of length na = ir(ai + a 2 )/2 with distance d between their axes. 
However, we would obtain a bad result if we simply multiply C by na in order to find 
Cmut- We should take into account the non-uniformity of the charge distribution along these 
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equivalent wires and then write C mut = CqP (where the effective length of two parallel 
charged wires P is smaller than na due to this non-uniformity). 

We use the King theory of loop antennas [16] in order to find effective length P. For it 
we write the Fourier expansion for the current induced in the reference loop by a magnetic 
field (with respect to the azimuth angle 0): 

J(0) = /(°) + /« cos + J {2) sin + . . . (35) 

It corresponds to the following expansion for the induced charges per unit angle: 

r(0) = Ti sin + r 2 cos + . . . (36) 

In the King theory it was proven that (for the case of small loop antenna {ka\ t 2 <C 1) with a 
finite impedance of the split) the omitted terms in expansions (35) and (36) are of the next 
order of smallness compared to the remained terms [16]. The uniform term represents the 
magneto-dipole mode of the induced current and has nothing to do with induced charges. 
Two other terms in (35) are responsible for the electric polarization of the loop and are 
related with the charge distribution. We can rewrite (36) as: 

r(0) = T max sin(0 + O ) (37) 



where T max = yr{ + r| . So, at it was already mentioned, half of the loop (0 O < < 7r + O ) 
contains positive charge Q, and the other half (— n + O < < 0o) contains negative charge 
— Q. The distribution of these charges is sinusoidal. The induced charges in the other loop 
have the opposite sign but the same sinusoidal distribution (see also Fig. 4). 

We obtain that the distribution of the charge along two effective equivalent parallel wires 
of length na is sinusoidal with zeros at the ends of these wires. It leads to the result P = 2a 
which gives together with (34) the approximate relation for the capacitance C mut : 

_ 2nee a 
arccosh J 

If the loop is excited by y— polarized electric field, the formula of the current distribution 
(35) keeps valid [16]. Coefficients jC - 1 ' 2 ) change in this case, but the charges +Q and — Q 
remain sinusoidally distributed on the two halves of each the loop. Therefore the formula 
(38) keeps valid, too. 
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Figure 1: Geometry of particles, (a) SRR particle (b) Omega particle 




(a) (b) 



Figure 2: Preparing the isotropic group of SRR particles, (a) A unit cell contains 6 particles. 
Particles located on opposite side of a cubic cell are rotated by 180°. (b) Composite can be 
prepared using the principle of chessboard. 
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(a) 



(b) 



Figure 3: Isotropic group of Omega particles (the same method as in 2). (a) Omega particles 
are positioned at the sides of a cubic unit cells, (b) The arrangement of unit cells shown 
here allows to obtain an effectively isotropic medium. 




Figure 4: Capacitive mutual influence of loops. The connection of equivalent capacitors C mut 
is shown. 
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Figure 5: Real part of the averaged magnetic polarizability of SRR versus frequency. The 
split width 5 is 0.1 mm. Analytical model (solid) and numerical model (dashed). 




4.5 5 
frequency (GHz) 



Figure 6: Imaginary part of the averaged magnetic polarizability of SRR versus frequency. 
The split width 5 is 0.1 mm. Analytical model (solid) and numerical one (dashed). 
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frequency (GHz) 

Figure 7: Real part of the averaged electric polarizability of SRR versus frequency. The split 
width is 0.1 mm. Analytical model (solid) and numerical model (dashed). 




frequency (GHz) 

Figure 8: Imaginary part of the averaged electric polarizability of SRR versus frequency. 
The split width is 0.1 mm. Analytical model (solid) and numerical one (dashed). 
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Figure 9: Real part of the averaged electric polarizability of SRR versus frequency. The split 
width is 0.4 mm. Analytical model (solid) and numerical one (dashed). 




Figure 10: Real and imaginary parts of the permittivity of the random SRR composite versus 
frequency. Analytical model (solid) and numerical model (dashed). 
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Figure 11: Real and imaginary parts of the permeability of the random SRR composite 
versus frequency. Analytical model (solid) and numerical one (dashed). 
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Figure 12: Real parts of the permeability (solid) and permittivity (dashed) of the SRR 
composite from cubic cells versus frequency. Analytical model. 
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Figure 13: Real parts of the effective permeability (solid) and permittivity (dashed) of an 
Omega composite formed by cubic cells versus frequency. The host permittivity is 2 and the 
length of a unit cell is 6.5 mm 
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Figure 14: Imaginary parts of the effective permeability (solid) and permittivity (dashed) of 
an Omega composite formed by cubic cells versus frequency. The host permittivity is 2 and 
and the length of a unit cells is 6.5 mm 
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Figure 15: Real part of the effective permeability (solid) and permittivity (dashed) of a 
random Omega composite. The host permittivity is 2. 
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